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Abstract 

Different methods are used to determine the scaling exponents associated with a time series 
describing a complex dynamical process, such as those observed in geophysical systems. Many of 
these methods are based on the numerical evaluation of the variance of a diffusion process whose 
step increments are generated by the data. An alternative method focuses on the direct evaluation 
of the scaling coefficient of the Shannon entropy of the same diffusion distribution. The combined 
use of these methods can efficiently distinguish between fractal Gaussian and Levy-walk time series 
and help to discern between alternative underling complex dynamics. 
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The evaluation of the scaling exponents is of fundamental importance to describe a num- 
ber of complex systems 0, 0, y|. The mathematical definition of scaling is as follows 0. 
The function $(r!,r 2 , . . .) is termed scaling invariant, if it fulfills the property: 

$(n, r 2 , • • •) = 7° $(7 6 n, 7 C ^2, • • •) • (!) 

Thus, if we scale all coordinates {r^} by means of an appropriate choice of the exponents 
a,b,c...., then we always recover the same function. This scaling invariance is the basic prop- 
erty that characterizes fractal functions 0,0,0]. The theoretical and experimental search 
for the correct scaling exponents is intimately related to the discovery of deviations from 
ordinary statistical mechanics. Fractal time series are particularly important in geophysics, 
as well as in several other field of research including biophysics and econophysics, where the 
phenomena of interest may present specific self-similarity patterns on different time scales. 

Two methods of analysis of time series commonly used to determine scaling properties 
are autocorrelation analysis and power spectral analysis [5]. The autocorrelation function 
of the fractal noise {^} results in the relation 



c ( r ) = < ^2 > ocr . (2) 



The power spectral representation of the same scaling property reads: 

/oo 
C(r) e- l27Tfr dr oc f 1 ' 2 " . (3) 
-oo 

It is easy to recognize the self-similarity or scaling property of the above two equations in 
their power-law form. The scaling exponent H was called the Hurst exponent by Mandelbrot 
in honor of the civil engineer Hurst who first understood the importance of scaling laws 
to describe the long-range memory in time series. In particular, Hurst was interested in 
evaluating the strength of the persistence of the annual level of the floods of the Nile river 
and, for such a scope, developed a time series analysis method to determine the scaling 
p_ H Q. 

A value H = 1 corresponds to 1/f-noise or pink noise. The adoption of a color name 
"pink" derives from the fact that a light source characterized by a 1/f spectrum looks pink. 
These type of noises are particular important because they represent a kind of perfect balance 
between randomness and order, or between unpredictability and predictability. In fact, for 
pink noises the autocorrelation function between two events separated by a time interval 



At = r is independent on r, (C(r) ~ const). Pink noises, if ~ 1, are found in countless 
natural phenomena from heart-beat intervals to music [7]. A value < H < 0.5 corresponds 
to antipersistent noise, H = 0.5 corresponding to uncorrelated or random noise, also known 
as white noise, and 0.5 < H < 1 corresponds to correlated or persistent noise. It is possible 
to extend the definition of H for values larger than 1. So, a value H = 1.5 corresponds to 
Brownian motion, which, as it is well known, describes the erratic motion of a particle, such 
as a pollen grain, in suspension on a fluid; this erratic motion is cause by random collisions 
between the particle and the molecules of the fluid $] . A value H = 2 corresponds to brown 
noise and a value H > 2 is known as black noise. These noises are characterized by a very 
smooth shape and may be adopted, for example, to generate artificial mountain landscapes 

0. 

It is important to point out that there are two common alternative deviations from 
ordinary statistical mechanics: anomalous Gaussian statistics and Levy statistics |9(. These 
two different statistics are indicative, in particular, of two different kind of complex noises: 

n 

the monofractal Gaussian noise and the Levy- walk noise 10J. These two types of noises 
present similar long-range correlation patterns, but are generated by quite different complex 
dynamics. The monofractal Gaussian noise, in its persistent form, presents long range 
memory in the sense that future events are strongly related to the frequency of occurrence 
of past events and the waiting time distribution between events has finite variance. The 
Levy-walk intermittent noise, instead, presents long-range correlation patterns which are 
generated by random waiting time intervals between events obeying to an inverse power 
law distribution with exponents that yield infinite variance, and there is no real correlation 
between events ll|. Figs. 1 show examples of these noises. 



Herein, we briefly describe two alternative time series scaling analysis methods, whose 
combined adoption can be used to distinguish the two above alternative noises. The diffusion 
entropy analysis (DEA) and standard deviation analysis (SDA) |l2|. Both techniques are 
based on the prescription that a time series {^} of N elements are the fluctuations of a 
diffusion trajectory jl2j ]. Note that there exist several other scaling analysis methods such 
as the detrended fluctuation analysis [l3( and several wavelet based methods 0,0, Q> 
which are variance based methods and are theoretically equivalent to the SDA. 



According to the prescription of Scafetta and Grigolini 



121, we shift our attention from 



the time series to the probability distribution function (pdf ) p(x, t) of the corresponding 
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FIG. 1: [A] Fractional Gaussian noise with H = 0.8; [B] Two forms of Levy-walk intermittent 
noise with ip(r) oc and \x = 2.5. B2 gives the frequency of impulses every 300 units of Bl. 

diffusion process, that is, the pdf of the diffusion process, p(x, t), is evaluated by means of 
the N — t sub-trajectories 

i 

X n {t)=J2&+" ( 4 ) 

8=0 

with n — 0, 1, . . . Therefore, x denotes the variable collecting the fluctuations and is referred 
to as the diffusion variable. The scaling property of the diffusion process, if it exists, takes 
the form 

where 5 is the scaling exponent. The DEA |l2| is based on the evaluation of the Shannon 
entropy S(t) using the pdf ((HJ). If the scaling condition of Eq. (JHJ) holds true, it is easy to 
prove that 

S(t) = - J p(x, t) ln[p(x, t)}dx = A + 5 ln(t) , (6) 

where, A is a constant. Numerically, we first evaluate the pdf with histogram of size-bin 
equal to the standard deviation of the data, and then use a discrete form of Eq. ( |B|) . 

The SDA 12j is based on the evaluation of the standard deviation D(t) of the same 
variable x, and yields 

D{t) = JT^^ (t)-(x-,t)} 2 a ^ (?) 

where (x; t) = J2n=o x n{t) is the mean value of {x n (t)}, and H is the Hurst exponent. 

If the data are fractal Gaussian noise the two scaling exponents are related to each other 
via the fractal Gaussian relation 

H — 5 . (8) 
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FIG. 2: DEA end SDA of: [A] a fractal Gaussian intermittent noise with ip( T ) ex P( — t /t) with 
7 = 25 and H = 5 = 0.75; the fractal Gaussian relation (|SJ) of equal exponents is fulfilled; [B] a 
Levy-walk intermittent noise with V ; ( T ) an d = 2.5; note the bifurcation between H = 0.75 

and 5 = 0.67 in accordance with the Levy- walk relation (|10l) . 

If the data are generated by a Levy-walk they are characterized by an inverse power law 
waiting time distribution of the type 

1 



ip(r) oc 



(9) 



[1 + tY 

where 2 < fi < 3, which ensures that although the first moment of r is finite, the second 
moment diverges. The scaling exponents are related to each other via the Levy-walk relation 
(LWR) Q 



0.5 < 5 
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12 



(10) 



There are several complex ways to generate a Levy- walk sequence, see Ref. [n, 
Some of these noises involve mixed Levy-Gaussian properties. The simplest Levy-walk 
sequence is a dichotomous signal made by a series of zeros and ones, where £ = 1 represents 
the occurrence of an event and £ = represents no event. The time intervals {rj} obeying 
to Eq. give the intervals between events. 

Thus, by evaluating 5 and H and using Eq. (JBJ and (fTTljl it is possible to distinguish 
the two kinds of time series while the adoption of only one of the two techniques can 
lead to a misinterpretation of the characteristics of a phenomenon. Figs. 2A and 2B show 
DEA and SDA applied to a fractal Gaussian intermittent noise with 5 = H = 0.75 and to 
a Levy walk intermittent noise with \x = 2.5, which correspond to 5 = 0.67 and H = 0.75, 



respectively 



HQ. 
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FIG. 3: [A] Waiting time distribution of flares. The distribution is fit with an inverse power law 
Eq. © with exponent fi = 2.12 ± 0.05 [Grigolini et al. 2002]. [B] SDA applied to the ACRIM 
composite TSI time series, sunspot number and cover sequence and global temperature anomalies. 
The uppermost line represents the theoretical Levy-walk scaling of the solar flare intermittency, 
H T = 0.94 ± 0.02, obtained via Eq. (JTUJ) with fx = 2.12 ± 0.05. The bottom line shows the scaling 
for random noise, H = 0.5, for comparison. 

Particularly interesting applications of the above scaling analysis techniques are found in 

i and global temperature patterns, 
Herein, we briefly summarize some 



geophysical phenomena such as earthquakes, solar flares and global temperature patterns, 



where long time series of data are available 
of our findings. 

In Ref. [TTJ it was shown that the waiting time interval distribution ip{ T ) between solar 
flares 2(| is an inverse power law of the type © with exponent fi = 2. 12 ±0.05. According to 
the LWR (|1U|) this would induce a Levy-walk with theoretical exponents H = 0.94±0.04 and 
5 = 0.89±0.04. Fig. 3A shows the waiting time distribution of flares. Fig. 3B shows the SDA 
applied to several solar data such as ACRIM TSI j^lj], sunspot cover , sunspot number 
[22j, global surface temperature anomalies (1856-2003) [24| and the theoretical prediction 
derived from the solar flare intermittency. The curves are quite parallel and suggest that a 
Levy like process regulate the dynamics of the solar activity and that the Earth climate seems 
to contain the same statistics. The latter statement seems further confirmed by Figs. 4A and 
4B that show DEA and SDA applied to the ACRIM TSI and the global temperature record 
showing the typical LWR bifurcation. We observe that if these findings are not accidental, 
they might imply the existence of a non-negligible complex Sun-Climate nonlinear coupling 
on a short time-scale as some studies seem to confirm 125, 2( 
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FIG. 4: [A] DEA and SDA applied to the ACRIM composite TSI time series. The two straight 
lines correspond to the scaling coefficients 5 = 0.88 ± 0.02 and H = 0.94 ± 0.02. [B] DEA and 
SDA applied to the global temperature anomalies (1856-2003) time series. The two straight lines 
correspond to the scaling coefficients 5 = 0.89 ± 0.02 and H = 0.95 ± 0.02. 



The above techniques can be applied also to earthquake concurrence [11]. An issue about 
seismic phenomena is whether: (1) they obeva statistics according to which the waiting 
times between Omori's earthquake clusters 12911 are uncorrelated from one another, as the 



traditional Generalized Poisson model 



27, m 



the "ETAS" model |30j; (2) the Omori's 
28| : (3) or whether the 



earthquake clusters obey to some knind of Levy-walk statistics 
data may also be characterized by intercluster 1/f long-range correlations between Omori's 
clusters that may disclose the earthquake conversations recently suggested by Stein j^ . 
Understanding the nature of the long-range correlation is fundamental for building reliable 
earthquake models. 

Fig. 5A shows the waiting time PDFs between earthquakes in California jsij using four 
earthquake magnitude thresholds M t = 1, 2, 3 and 4. The PDFs show an initial Omori's law 
2^| (.P(t) oc 1/t), but the pdf tails present a large inverse power law exponent [i > 4 and 
may even approach an exponential (or Poisson) distribution asymptotically. The Omori's 
law is determined by the short-range correlated aftershocks |2i| and lasts for a time that 
increases with the magnitude threshold. Fig. 5B shows the DEA and SDA applied to the 
intermittent time signal £(£), where t is the physical time, obtained by assigning a value 
equal to 1 at the occurrence of an event, and a value equal to when no event occurred. 
The latter figure suggests that the data fulfill FGR ©. Thus, beyond the Omori's law, the 
earthquake clusters might be uncorrelated if the observed super-diffusion 5 = H = 0.94 is 
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FIG. 5: [A] pdf of the waiting times Tj of earthquakes with a magnitude M > Mt = 1, 2, 3 and 4. 



The initial P{t) oc 1/t is the Omori's law |29|. [B] DEA and SDA of the intermittent time signal 
for magnitude M > Mt = 1. The data are fitted with scaling exponents 5 = 0.944 ±0.008 and 
H = 0.943 ± 0.004. The uppermost solid line with H = 0.97 corresponds to the expectation of H 
if the Levy-walk condition 1)10(1 holds true. 

generated by a long-tailed Omori's law involving multiple clusters, or there might be the 
possibility that the clusters are correlated as a 1/f Gaussian noise 11 ill. In the latter case, 

Tin 

traditional earthquake models such as the Generalized Poisson model 27, 28 or the "ETAS" 

n 

model [3fJ should be improved by adding additional correlations between clusters. 

In conclusion, we have discussed some properties of complex time series analysis showing 
two different types of anomalous statistics: fractal Gaussian noise and Levy-walk noise. 
We have shown how the multiscaling comparative analysis of time series can be used 
to distinguish the two types of noises and applied it to study some complex patterns of 
geophysical phenomena. Thus, we conclude that there are some difficulties in interpreting 
intermittent sequences. Models with alternative statistics can reproduce some patterns of a 
time series equally well. This fact suggests the need of an analysis involving complementary 
tests for addressing complex systems. 



Acknowledgment: N.S. thanks the ARO for support under grant DAAG5598D0002. 



[1] B.B. Mandelbrot, The Fractal Geometry of Nature, Freeman, New York, (1983). 
[2] Feders J., Fractals, Plenum Publishers, New York, (1988). 



S 



[3] Peitgen H.-O., J. Hartmut, D. Saupe, Chaos and Fractals, new frontiers of science, sec. edition, 

Springer, New York, (2004). 
[4] Goldenfeld N., Lectures on Phase Transitions and the Renormalization Group ( Perseus Book, 

Reading, Massachusetts, 1985). 
[5] R. Badii and A. Politi, Complexity, Hierarchical structures and scaling in physics, Cambrige 

University Press, UK, (1997). 
[6] Hurst H.E., R. P. Black, Y.M. Simaika, LongTerm Storage: An Experimental Study, Constable, 

London, (1965). 

[7] Schroeder M., Fractals, Chaos, Power Laws: Minutes from an Infinite Paradise, W.H. Freeman 
& Company,(1992). 

[8] Gardiner C.W., Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sci- 
ences, 2 nd edition, Springer- Verlag, New York, New York, (1997). 

[9] A.I. Khinchin, Mathematical Foundations of Statistical Mechanics, Dover Publications, Inc. 
New York, (1949). 

[10] Shlesinger M.F., B.J. West, J. Klafter, "Levy dynamics of enhanced diffusion: Application to 

turbulence," Phys. Rev. Lett. 58, 1100 (1987). 
[11] Scafetta N., and B.J. West, "Multi-scaling comparative analysis of time series and a discussion 

on earthquake conversations in California," Phys. Rev. Lett. 92 138501 (2004). 
[12] Scafetta N., and P. Grigolini, "Scaling detection in time series: diffusion entropy analysis," 

Phys. Rev. E 66, 036130, (2002). 
[13] Peng C.-K., S.V. Buldyrev, S. Havlin, M. Simons, H.E. Stanley, and A.L. Goldberger, "Mosaic 

organization of DNA nucleotides," Phys. Rev, E 49, 1685, (1994). 
[14] Mallat S.G., A Wavelet Tour of Signal Processing (2nd edition), Academic Press, Cambridge 

(1999). 

[15] Muzy J.F., E. Bacry, A. Arneodo, "The Multifractal Formalism Revisited with Wavelets," 
Int. J. Bifurc. Chaos 4, No. 2, 245-302 (1994). 

[16] Percival D.B., and A.T. Walden, Wavelet Methods for Time Series Analysis, Cambrige Uni- 
versity Press, Cambrige (2000). 

[17] Grigolini P., D. Leddon, N. Scafetta, The Diffusion entropy and waiting time statistics of hard 
x-ray solar flares, Phys. Rev. E 65, 046203 (2002). 

[18] Scafetta N., and B.J. West, Solar Flare Intermittency and the Earth's Temperature Anomalies, 



9 



Phys. Rev. Lett. 90, 248701 (2003). 
[19] Scafetta N., V. Latora and P. Grigolini, "Scaling without detrending: the diffusion entropy 

method applied to the DNA sequences," Phys. Rev. E 66, 031906 (2002). 
[20] Solar Flare Catalog, (2003). http://umbra.nascom.nasa.gov 

[21] Willson R.C., and A.V. Mordvinov , "Secular total solar irradiance trend during so- 
lar cycles 21-23," Geophys. Res. Lett., 30, 1199, doi: 10.1029/2002GL016038 (2003). 
http://www. acrim. com 



[22 
[23 

[24 
[25' 

[26 



[27 
[28 

[29 



[30 
[31 



Solar Influences Data analysis Center (2003). http://sidc.oma.be/index.php3 



Royal Greenwich Observatory/USAF/NOAA Sunspot Record 1874-2003 (2003). 

http://science.msfc.nasa.gov/ssl/pad/solar/greenwch.htm 

Climatic Research Unit, UK, (2003) http://www.cru.uea.ac.uk, 

D'Aleo Joe, "Solar influence may be evident in three of the last four seasons," (2002). 
http:/ '/www.intellicast. com/DrDewpoint/Library/1339 

Jackman C.H, and R.D. McPeters, "The effect of solar proton events on ozone and other 
constituents," Chapter in Solar Variability and its Effects on Climate by Pap J.M, and P. 
Fox, Geophysical Monograph Series Volume 141 (2004). 

Stein R.S., "Earthquake Conversation," Scientific American, 288 Jan 72 (2003). 
Mega M.S., P. Allegrini, P. Grigolini, V. Latora, L. Palatella, A. Rapisarda and S. Vinciguerra, 
"Power-Law Time Distribution of Large Earthuakes," Phys. Rev. Lett. 90 188501 (2003). 
Omori F., "On the aftershocks of earthquakes," J. College Sci. Imper. Univ. Tokyo 7, 111 
(1895). P. Bak, K. Christensen, L. Danon and T. Scanlon, "Unified Scaling Law for Earth- 
quakes," Phys. Rev. Lett. 88 178501 (2002). 

Kagan Y.Y., and L. Knopoff, "Statistical short-term earthquake prediction," Science 236, 
1563 (1987). 

Southern California Earthquake Data Center (2003). http://www.scecdc.scec.or g/ftp/cata logs/SCNS/ 



10 



